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We study the dynamical evolution of a phase interface or bubble in the context of a \(f> 4 + g(f> G 
scalar quantum field theory. We use a self-consistent mean-field approximation derived from a 2PI 
effective action to construct an initial value problem for the expectation value of the quantum field 
and two-point function. We solve the equations of motion numerically in (l+l)-dimensions and 
compare the results to the purely classical evolution. We find that the quantum fluctuations dress 
the classical profile, affecting both the early time expansion of the bubble and the behavior upon 
collision with a neighboring interface. 
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I. INTRODUCTION 

During first order phase transitions, bubbles or do- 
mains of the lower free energy phase (true vacuum) are 
nucleated in a metastable, or false vacuum, phase. Even 
at zero temperature, bubbles are induced by quantum 
effects, but they may also be thermally activated. The 
theory of droplet formation, describing the onset of nu- 
cleation, is by now venerably old and well established. 
It stretches back to the work of Becker and Doring 
and Langer Q in statistical physics and later to that of 
Coleman || in the context of relativistic quantum fields, 
among others |Q. 

The theory of droplet nucleation, although successful, 
leaves almost all dynamical questions unanswered: what 
happens to the system once the bubble is nucleated? The 
general phenomenology of bubble expansion and coales- 
cence must address how a semi-classical field solution (the 
bubble) propagates in the presence of quantum or ther- 
mal fluctuations for long times, i.e. how these fluctua- 
tions interact quantum mechanically with the interface, 
and how the full self-consistent system may be described 
classically by hydrodynamics, for example of front prop- 
agation in media. 

All of these questions can be easily posed and are, in 
principle, answerable in the context of quantum field the- 
ory. Tackling them quantitatively however requires a 
combination of non-perturbative analytical and numer- 
ical techniques that are just now beginning to emerge. 
The aim of the present paper is to take the first steps to- 
ward studying the nucleation and dynamical propagation 
of bubbles together with their self-consistent quantum 
fluctuations in relativistic quantum field theory. 

The theory of droplet nucleation tells us that there 
are subcritical bubbles which decay away and also super- 
critical bubbles which feed on the energy released by the 
phase transition to grow until the true vacuum phase has 
obliterated the false vacuum entirely. Coleman dubbed 
this the fate of the false vacuum pi . 



Details of the dynamics described heuristically above 
are notably absent. To address the question of what is 
the critical bubble size, what is the shape or profile of the 
bubble as it expands, and whether the bubble wall experi- 
ences viscous drag, we need a thorough understanding of 
the nonequilibrium quantum field dynamics. A classical 
analysis based on global properties of Lorentz invariance 
is present in Q , but this of course leaves out fluctuations 
and hence both virtual or real particles. 

The inclusion of (self-consistent) particles or fluctua- 
tions leads to a panoply of new phenomena that must 
be considered for the complete description of the phase 
transition. Recognizing this, Coleman left a number of 
open questions concerning the effects of fluctuations (par- 
ticles) on interfaces and vice- versa ||. These issues are 
not manifest in the strict context of the nucleation prob- 
lem. 

The first of Coleman's questions is, what happens when 
a bubble encounters particles? This phenomenon is cen- 
tral to scenarios of early Universe baryogenesis and has 
been addressed in this context to some extent || . Baryo- 
genesis remains the most important motivation for the 
study of bubble wall dynamics in relativistic settings. 
Several works have recently addressed the problem 
of computing the asymptotic velocity and shape of Higgs 
field bubbles at temperatures near the electroweak phase 
transition. These approaches treat the bubble wall as a 
classical field background immersed in a bath of thermal 
fluctuations which obey effective transport equations for 
their occupation number distributions. This treatment is 
appropriate if the bubble wall moves sufficiently slowly, 
contains only "soft gradients" , and if quantum coherence 
is unimportant. Thus a transport approach will neces- 
sarily fail at sufficiently low temperatures and/or under 
severe supercooling. In these more difficult cases, the 
direct field theoretical methods developed here become 
essential. Quantum first order phase transitions in non- 
relativistic systems Q may provide an interesting lab- 
oratory for testing the non-relativistic analogue of the 
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zero-temperature methods described below. 

Coleman's other questions are concerned with the pos- 
sibility that bubbles may be induced by fluctuations (and 
perhaps even created at particle scattering experiments 
H ) and with particle production resulting from the colli- 
sion of two bubble walls. Both phenomena necessitate a 
dynamical non-perturbative treatment of quantum field 
theory valid for long times. For this reason they have 
remained poorly understood. 

In recent years, the availability of numerical methods 
to solve for the time evolution of quantum fields has given 
rise to a resurgence of interest in such problems. The 
causal formalism suited to initial value formulations of 
field theory dynamics has been employed in various ap- 
proximation schemes in an effort to isolate the relevant 
features of a quantum kinetic theory from first princi- 
ples§J|4|. 

In this paper we consider a scalar quantum field theory 
which exhibits a first-order phase transition. Assuming 
that the field is in the "false vacuum" before the transi- 
tion and is brought out of equilibrium by the nucleation 
of bubbles in this "true vacuum" phase, we study the 
detailed dynamics of bubbles which we impose as initial 
conditions. Because of the computational effort required 
in the quantum theory, we restrict our attention to (1+1)- 
dimensional spacetime. 

We consider the purely classical field evolution as well 
as a self-consistent quantum evolution in the Hartree 
approximation at zero temperature. The generalization 
of the formalism to include both higher-order interac- 
tions |l5| ] and/or finite temperature is straightforward. 
Thermal effects lead to qualitatively different physics, 
and we intend to analyze these physical consequences in 
detail in a future work. 

The main results of this paper are: 1) Whether a true 
vacuum bubble is critical is determined by the extrem- 
ization of the energy, not the action. While this should 
be obvious from the point of view of an initial value prob- 
lem, there has traditionally been some confusion of the 
critical (spacetime) radius for the bounce, Rb with the 
critical (purely spatial) radius for growth, which we label 
Re- The two values are related by a constant of propor- 
tionality 

Re = (d-1) 
R B d 

in d spatial dimensions; hence for (l+l)-dimcnsions, any 
bubble is critical. More precisely, the critical bubble size 
is constrained in one dimension only by the thickness of 
the bubble. 2) The bounce determines the correct profile 
of the bubble wall, but induced, super-critical bubbles 
with larger or smaller radii still grow and asymptote to 
shifted light-cones. The bounce solution is unique and 
identifiable in that it asymptotes to the light cone from 
the origin. These results are already manifest in the clas- 
sical description. 3) Including quantum effects at the 



level of the Hartree approximation does not change the 
qualitative features-constrained by Lorentz invariance-of 
bubble growth at zero-temperature. However, in much 
the same sense that quantum fluctuations render the 
quantum effective potential different from the classical 
one, they do affect the detailed shape of the bounce. 
Hence, 4) The proper description of quantum bubble 
dynamics necessitates a self-consistent bounce which in- 
cludes a prescription of the quantum fluctuations at the 
time of nucleation. 5) The behavior of colliding bub- 
bles does indicate a qualitative difference between the 
classical and quantum behavior. In our model, the clas- 
sical bounce appears remarkably stable against bubble 
coalescence-exhibiting elastic collisions off neighboring 
expanding bubbles for very long times (possibly forever). 
The quantum evolution, on the other hand, displays the 
more expected behavior wherein the bubbles disappear 
by transferring energy to intermediate frequencies on a 
time scale of the same order of magnitude as the bubble 
size. 

In Section II we summarize the semiclassical theory 
and its predictions; we present our model and highlight 
some of the classical dynamical details which emerge 
anew due to the consideration of the dynamics as an 
initial value problem and due to our specialization to 
(l+l)-dimensions. In Section III, we extend the analysis 
of the dynamics to a self-consistent Hartree-like approx- 
imation and discuss the simplifications it involves. We 
summarize our results for the propagation and collision 
of bubble walls in the quantum theory in Section IV. 
In Section V we discuss the many interesting possibili- 
ties for the application of the methods of this paper to 
related questions as well as the refinements necessary to 
render the long-time evolution of self-consistent quantum 
fluctuations more realistic. 

II. (SEMI) CLASSICAL DYNAMICS 

Ultimately, we would like to understand the quantum 
dynamical evolution of a generic bubble of true vacuum 
(induced perhaps by coupling to other fields or sources). 
We should naturally do first what we can in the classical 
regime, where we may apply the literature on semiclassi- 
cal field theory methods in the bubble nucleation prob- 
lem and connect to other numerical studies |l(]] . The 
relativistic picture was elegantly framed by Coleman in 
Ref. ||, so we shall parallel that analysis, working out a 
specific example in full detail. 

Coleman set out to compute the decay rate of the false 
vacuum in a scalar theory described by the Lagrangian 

£ = \(d,cf>) 2 -V(rt . (1) 

By analogy with the semiclassical analysis of barrier pen- 
etration, he obtained the exponent in the vacuum decay 
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rate in terms of an instanton solution of Euclidean space- 
time which he called the bounce. The bounce function is 
a saddle point of the Euclidean action, 



S E [<t>] 



d d x dr 




hence it satisfies Euclidean "equations of motion". The 
solution is subject to appropriate boundary conditions at 
the origin and at infinity. It is understood that such a 
function will always exist and will be 0(D) invariant for 
D = d+ 1 spacetime dimensions, thus depending only on 
a radial coordinate p = \/t 2 + x 2 , 0& = 4>b(p)- Then the 
bounce equation and corresponding boundary conditions 
take the form 



2 </> b D - 1 d<t> b _ SV(<t> b ) 



dp 2 



dp 



lim <f>(p) = fa 

p — >oo 



and 



dtp 
dp 



(3) 
(4) 



p=0 



A happy consequence of the bounce solution is that the 
real-time classical equation of motion is just the analytic 
continuation of the bounce equation to real time and thus 
is solved by the analytic continuation of the bounce from 
Euclidean spacetime to Minkowski, i.e. with r — * it, 



fa(p) -> MVx 2 - 1 2 ). 



(•5) 



The shape of the bounce becomes the profile of the bub- 
ble wall, and this shape remains unaltered as the wall 
describes a hyperbola in spacetime, reaching the speed 
of light asymptotically. Lorentz invariance allows us to 
solve for the bubble profile in all space and time and not 
as an initial value problem. 

A number of features of the spontaneous decay of the 
false vacuum are severely constrained. The most favor- 
able shape of the bubble is constrained by stationariza- 
tion of the action. That its asymptotic velocity is c (or 
1 in natural units) is enforced by 0(D) invariance of the 
bounce in Euclidean spacetime or 0(D — 1, 1) Lorentz in- 
variance. Thus it is impossible in the absence of Lorentz 
violating fluctuations for the bubble wall to experience 
drag. Classically, any drag would have to be put in by 
hand in the equations of motion^. 

Let us make the analysis concrete by specializing to a 
potential of the form 



*We have verified that the addition of a simple drag term 
rjdt(p in the dynamical field equations does indeed result in an 
asymptotic interface velocity smaller that that of light. The 
study of the stochastic bubble wall trajectory in the presence 
of phenomenological damping and noise is in itself an inter- 
esting problem, see Jl7| . 



V[ct>] = -m 2 2 + A0 4 + . 



(6) 



+ v(fa\, (2) 



For the purposes of comparison with an analytic, thin- 
wall approximation, it is convenient to have a potential 
with degenerate minima in a simple parametric limit. We 
can rewrite the potential in terms of the parameters (/>q 
and 7 such that 



V[4>] =9^ (4? - 4>lf - i<t> 2 



(7) 



where (b 2 



-A/25 and 7 = (A 2 /.g - 2m 2 ) /A. The thin 
wall approximation relies on neglecting the term pro- 
portional to 7. In practice the thin wall obtains in the 
limit where the energy difference between the two min- 
ima e = V(fa-) — V(fa.) — 7</>q vanishes. As a working 
example we use the parameters m 2 — 4, A = —0.8 and 
g = 0.07. Then 7 = 0.29. With these choices the poten- 
tial for the bounce calculation is shown in Fig. [j]. Also 
shown is the degenerate potential (7 = 0). 



V{<t>) 




FIG. f . The classical potential (lower curve) for the bounce 
calculation with minima at <p- = and </>+ = 2.46. The 
upper curve is the degenerate potential (see Eq. (^)), used to 
compute in the thin wall approximation. Its second minimum 
occurs at d> = 2.39. 



We construct the bounce solution in the thin wall ap- 
proximation, following Ref. [^), and by direct numerical 
integration. The thin wall bounce consists of solving for 
the configuration of </>, piecewise in p such that 



b- = <po + 0(7), 

hmu(p-R), P 



P< R 
t R 



(8) 



= 6+ = 0, 



p > R. 



Here </> wa ii is computed in the degenerate potential, i.e. 
by solving 



d? </>wall 

dp 2 



= U'[4> 



'wall J 



U[4>] =g^[^ + — 



(9) 
(10) 
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The solution has been obtained previously pi 



^wall(P) 



1 + ew ' 



(11) 



where [i = y/Sgcf) 2 , — 4.21, and p? is the second deriva- 
tive of the potential U"[<fio] evaluated at </>o. Up to a 
correction of order e, it is the mass of excitations around 
the true minimum. To complete the thin wall approxi- 
mation R is determined variationally, as the value that 
extremizes the Euclidean action. In the thin wall ap- 
proximation to the action, the interface itself amounts to 
a surface term, whereas the contribution from the two 
piecewise constant parts is proportional to the bubble's 
spacetime volume. Hence in (l+l)-dimensions 



Se 



-irR 2 e + 2nRa, 



(12) 



where a is a surface tension. The stationary point of (|lz 
is equivalently the zero-energy value of R in configuration 
space: 



Re + a = 



so that R = a/e with 
dx 



1 / #wall 

2 V dx 



U[</> 



wall 



(13) 



(14) 



With our parameters, a — 3.05, and R = 1.82. The two 
solutions, the exact bounce computed numerically and in 
the thin wall approximation, are shown in Fig. g. 



FIG. 2. The numerical solution for the bounce profile 
(points) and the thin wall approximation from Ref. jl8| . The 
thin wall approximation agrees well with the exact numerical 
computation where it should, i.e. at p ~ R c — 1.82. 

We note that the thin wall bounce is a good approxi- 
mation to the exact solution in the neighborhood of the 
inflection point i.e. at p ~ R c — 1.82. The true con- 
dition for the validity of the approximation H is that 
piR 3> 1. With our parameters, p = ^/8g<j> 2 = 4.21, so 
that pR = 7.66 ^> 1, which is shy of an order of magni- 
tude larger than unity. In fact, as depicted in Fig. 0, the 
bubble wall does not appear very "thin" . 



To confirm the predictions of the Euclidean solution 
and test our numerical methods we now solve the clas- 
sical real time equations of motion (i.e. in Minkowksi 
space) for the field <p explicitly. As initial conditions, we 
use a bubble at rest with the profile given by the approx- 
imate analytic form ([ll]). By following one point on the 
bubble wall with time, we observe the predicted trajec- 
tory, shown in Fig. |3[ 



FIG. 3. Spacetime trajectory for a (l+l)-dimensional bub- 
ble started at rest at t = is well described by the hyperbola 
x = VR 2 + t 2 with R = 1.82. 

We have so far said nothing about induced vacuum 
decay, in which a bubble may be nucleated with an ar- 
bitrary size and shape. As seen in Eq. ([l3|), the sponta- 
neous decay of the vacuum via the bounce costs exactly 
no energy. However, if the energy of an initial bubble 
configuration may be accounted for from another source, 
it need not be zero. We would like to know how such a 
bubble will evolve given what we know about the bounce. 
We first consider varying the size of the bubble without 
changing its interface profile. 

Generally the fate of an induced bubble-growth or 
decay-can be determined by considering the energetics 
of the corresponding initial value problem. In this case, 
we do not have a solution for all spacetime, but only a 
spatial profile at some initial time t = 0. We will as- 
sume that the bubble wall is initially at rest, i.e. that its 
kinetic energy is zero f\ The total energy is given by 



E = / d d x \ (d t <(> h ) 2 



i(V0 b ) 2 + U(0 b ), (15) 



where the static part can again be computed in the thin 
wall limit. It is Static = a Sd ~ where Sd and Vd 
are the surface and volume of the bubble in d = D — 1 
spatial dimensions. The energy has a local maximum at 



^This is a weak assumption from the point of view of the 
initial value problem. From the Euclidean spacetime analysis, 
it is a consequence, not an assumption 
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Re = (d — l)c/e (recall that the bounce has an action 
extremum at Rg = dcr/e). For the interface to move 
while globally conserving energy it is necessary that it can 
lower its static energy, so that the difference is converted 
into kinetic energy. Bubbles with R > Re do so by 
growing while those with R < Re must shrink. Hence 
Re defines the critical radius for growth. A particular 
case is that of the bounce which corresponds to the choice 
of Rb which enforces Static = 0. The bounce always 
falls in the class R > Re and therefore always grows. 
However, it is not the case, as is frequently assumed in 
the literature, that the bounce radius is the critical radius 
as defined by the onset of dynamical growth. We plot an 
example of the energy in 3 spatial dimensions in Fig. ^(a). 




(b) d = 1 

FIG. 4. The dependence of energy and action on the radial 
coordinate in the thin wall approximation: (a) a "spherical" 
bubble in 3d. The critical radius for growth Re is the maxi- 
mum of the energy and is smaller by a factor of (d — 1) jd than 
the radius Rb which maximizes that action and corresponds 
to the zero of energy, (b) In Id, the critical radius Re goes 
to zero. 

The case of (l+l)-dimensions (d = 1 above) is special 
because the interface energy saturates: once the bubble is 
larger than twice the thickness of the wall, the "surface" 
term no longer depends on the size of the bubble. Since 
the bubble can gain kinetic energy by converting false 
vacuum, a thin wall bubble will always grow in (1+1)- 
dimensions, i.e. Re = 0. This is seen clearly in Fig. ||(b), 
where the thin wall approximation is again assumed. If /i 
is a measure of the inverse thickness of the bubble, then 
the thin wall approximation amounts to pR 3> 1. In 
practice, the finite size of the wall profile constrains fi so 
that R c ~ 1/fi gives a lower limit on the critical bubble 
size. 

We verified this behavior numerically as well as the 
fact that the thin wall prediction of R E =2 — Rb =1 = 



1.817 is an excellent approximation to the real critical 
value for bubble growth in d = 2, which occurs between 
R = 1.810 — 1.815. Changing the radius amounts to using 
4>{p — R) with the analytic form ( |lT| ) and variable R as an 
initial condition. We show the trajectories obtained for 
bubbles with different initial sizes in Fig. |[ Note that the 
trajectories for bubbles with radii smaller or larger than 
the bounce asymptote to shifted light cones. Naively, one 
might have expected all solutions to approach the light 
cone from the origin. 
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FIG. 5. Spacetime trajectories for the motion of several 
(l+l)-dimensional interfaces started at rest at t = and 
with different bubble radii. All trajectories are well described 
by hyperbolae x = ±(xo + \/R 2 + t 2 ), with origins set at 
t = 0. Thus the interface velocity at large times approaches 
the speed of light as in the case of the bounce, where xo = 0. 

In order to explain these trajectories, we explicitly con- 
struct a piecewise solution for all (real) time with an arbi- 
trarily initial size. The real time equations of motion are 
invariant under a coordinate shift x — > x + Xq. Consider 
cutting the function <j>b(x,i) along x — at t — and 
shifting the positive x piece by a positive amount xq and 
the negative x piece by —xq. In the gap created in the 
center, extend the value of <f>(x = 0) to be constant. This 
function looks exactly like the bounce at the interfaces 
but has more true vacuum sandwiched in the middle, 
see Fig. |6[ The boundary conditions on the bounce (Q) 
ensure continuity up to and including first derivates at 
±Xo- (Excising a homogeneous region in the center, one 
can similarly generate a smaller bubble.) 



t>0 













static vacuum 




, t>0 








-X 







FIG. 6. Schematic construction at t = of a solution to 
the real time equations of motion for all i > 0. 
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A constant field in true vacuum is a static solution 
of the equations of motion, hence the piecewise solu- 
tion will be static in the region |a;| < xq and expand 
like the bounce outside this region. Moreover, we recall 
that although R is determined variationally in the thin 
wall approximation-this freedom is exploited by our ini- 
tial condition-solving the bounce equation directly (un- 
der the assumption of O(D) symmetry) determines R 
uniquely. The coordinate shift transforms 



4>{^x 2 - t 2 - R) -» <j>(y/(x T x ) 2 -t 2 -R) , (16) 

so the trajectory of the interface is given by x = ±(xo + 
\JR? + 1 2 ) in real time, implying that bubbles induced 
with different sizes will asymptote to different light cones 
as in Fig. ||. 

So far we have altered the bubble size but not its pro- 
file. Since to distort the shape would be to perform a 
variation in function space, it is much more difficult to 
quantify the difference in the evolution of an arbitrary 
bubble from that of the bounce. Hcuristically, we under- 
stand that if a bubble starts out with the wrong profile, 
it will try to change its shape to the bounce before it 
expands as described above. The critical size argument 
can be affected by this: since a small bubble shrinks as 
it deforms, a slightly supercritical bubble with the wrong 
shape may still collapse away. A large initial bubble on 
the other hand will convert the false vacuum energy into 
kinetic modes which both expand and distort the bubble. 
We have observed both types of behavior in numerical 
evolutions. 



III. QUANTUM BUBBLES 

Extending the analysis of the previous section to quan- 
tum fields introduces a host of new challenges. From the 
point of view of dynamical equations of motion, we must 
now consider operator equations, which in fact translate 
into an infinite Dyson-Schwinger hierarchy of equations 
for the n-point Green's functions of the theory. For prac- 
tical purposes, the hierarchy must be truncated somehow, 
and this truncation introduces an approximation, often 
in the form of a self-consistent ansatz |l9| . Furthermore, 
even upon keeping a finite number of connected correla- 
tion functions the resulting system of coupled equations 
still describes an infinite number of degrees of freedom 
with non-linear interactions. Such a system must in gen- 
eral be solved numerically in a finite computer, which can 
lead to artifactual effects. Thus it is of utmost impor- 
tance to identify practical approximation schemes that 
capture the qualitative essence of the dynamical evolu- 
tion of a quantum field theory as compared with its clas- 
sical counterpart. This challenge is not devoid of uncer- 
tainty since not much is known about the evolution of 
truly quantum many body systems far away from ther- 
mal equilibrium. 



The formulation of the problem can be cast into one 
formalism which has been developed and explored in sev- 
eral works and has resurfaced in recent years with re- 
newed vigor |ll],[l^,[l4| . It is based on a two-particle ir- 
reducible (2PI) effective action for the field and the two- 
point function T[4>, G] formalized by Cornwall, Jackiw 
and Tomboulis (CJT) @. The Schwinger-Keldysh 
closed time path (CTP) is employed to make causality 
explicit through an appropriate real time prescription in 
the path integral time contour and associated Green's 
functions; hence it is frequently referred to as the CJT or 
the 2PI-CTP formalism. It is in the evaluation of T[(f>, G] 
that some sort of expansion (in loops or in powers of the 
coupling constant or in powers of 1/N for example) is 
carried out. 

In this paper we shall employ the 2PI-CTP formal- 
ism at the level of the Hartree approximation, a mean- 
field approximation which amounts to keeping the 2PI 
diagrams which are lowest order in coupling constants 
(bubble diagrams). The Hartree approximation is well 
understood to be equivalent to a Gaussian variational 
ansatz in the Schrodinger functional formalism p0[ and 
results in Hamiltonian dynamics pjj . It is also qualita- 
tively similar to the systematic large- A approximation at 
leading order in 1/N. For our purposes, A = 1 and the 
Hartree approximation is motivated chiefly because it is 
extremely simplifying. A large N scalar theory is also 
inappropriate for the description of first order transition 
dynamics as it results invariably in second order criti- 
cal phenomena (unless the O(N) symmetry is explicitly 
broken, whence criticality is erased and the transition 
becomes an analytical crossover). 

Time evolution of spatially inhomogeneous quantum 
fields, even at this level of approximation has only re- 
cently been produced numerically (2^-|^] . To our knowl- 
edge, the growth of a critical self-consistent quantum field 
bubble has not been previously demonstrated. 

The starting point for the 2PI-CTP formalism is the 
effective action T[<f>, G] of 0: 

r[cj>,G] = S[4>] + ^TrlnG" 1 + -iMT^G 

+r 2 (0,G) + const. , (17) 

where iD^ 1 (4>) is the classical inverse propagator 



iD 1 (x,y;<f>) = 



6S[<f>] 



8(t>{x)8(t>{y) 



(18) 



and T2 sums the 2PI vacuum to vacuum diagrams with 
propagators set to G and interaction vertices obtained 
from the shifted Lagrangian £[</>—> <fi + <p]. 



G 




FIG. 7. Vacuum bubble diagrams in the Hartree approxi- 
mation for the 2PI effective action of our model. Lines denote 
the full propagator G in the nontrivial (f> background. 

In the Hartree approximation, we consider 2PI dia- 
grams in Fig. [?], hence 



T 2 = -3(A + 15 5 </> 2 ) / d 2 x G(x, xf 

-15.9 / d 2 xG{x,x) 3 



(19) 



We observe that the two-point function enters into the 
Hartree approximation only as a local function evaluated 
at one spacetime point. This is apparent from the di- 
agrams in Fig. Q which, upon opening a line, become 
proper self-energy diagrams. The inclusion of any dia- 
gram with two or more vertices would introduce non-local 
kernels in the equations of motion. 

The stationarity conditions on the effective action in 
the absence of sources, 



Sr&G] _ 6T[<j>,G\ 
86 5G 



= 0, 



(20) 



lead to equations of motion for the field 



(□ + m 2 + AX6 2 + 6g6 4 + 12(A + 5# 2 )G(x, x) 

+90.g / d 2 x G(x, a;) 2 ) = 0, (21) 



and for the Wightman two-point function 

(o + m 2 + 12X6 2 + 30g6 4 + 12(A + 15# 2 )G(x, x) 

+90gG{x,x) 2 ^JG(x,x') = . (22) 

Let us remark that there are logarithmic and quadratic 
divergences in the bare theory in the self-energy and 
energy-momentum tensor respectively. However, any 
scalar field theory in (l+l)-dimcnsions is renormalizablc 
using the standard countcrterm procedure. For the 6 6 
interaction, this has been carried out explicitly in an 
early analysis of the quantum effective potential P5[ : the 
logarithmic divergence from the scalar loop in (1+1)- 
dimensions affects both the mass and the quartic cou- 
pling constant A, which need to be renormalized. In 
practice, even though we solve the equations of motion 
numerically on a discretized spatial lattice, we introduce 



time-independent counterterm subtractions as renormal- 
ization conditions for the "physical" mass, quartic cou- 
pling and energy density in vacuum. We fix the physical 
(dressed) mass to coincide with the bare mass of the clas- 
sical theory, i.e. we employ a subtraction (detailed be- 
low) which sets the two-point function to zero in the false 
vacuum. The quadratic divergence in the energy density 
is just the usual zero-point energy of vacuum, and again 
we choose to set the energy density in the false vacuum 
to zero. 

In order to solve the system of equations, we need to 
specify initial data for the field and the two-point func- 
tion. This is more or less straightforward for the field 
expectation value, which we identify with the classical 
field bubble profiles of Section ||. Thus, as an ansatz, 
we initialize the quantum field expectation to have the 
same shape as the bounce function described in Eq. jTl]). 
The initial specification of the two-point function, how- 
ever, is informed by the full details of the noncquilibrium 
ensemble. As a starting point, we use the most naive 
initial conditions possible, which is to ignore the forma- 
tion of the initial bubble completely. In other words, 
the two point function is initialized for a zero- or finite- 
temperature distribution about the false vacuum every- 
where in space, including inside the bubble. We are in 
the process of improving this prescription and comment 
on this below. 

We proceed by decomposing the equal-spacetime prop- 
agator in a mode basis. From here on we use x to de- 
note the space coordinate alone instead of spacetime, i.e. 
x^ = (x,t) and G(x ti ,x^ 1 ) can be written as 

G(x,t) =^|<Mx,i)| 2 (2n s (fc) + l) , 

k 

where ns{k) is the Bose-Einstein distribution 

1 



(23) 



n B = 



1 



(24) 



The mode functions ipk are initialized in a plane wave ba- 
sis with the mass of excitations around the false vacuum: 



tpk(x) 



Ak-x 



with 



Vk 2 



m 



(25) 



(26) 



The mass in the dispersion relation is self-consistently 
dependent on the values of the mean field and fluctua- 
tions. The log-divergent, zero-temperature part of G is 
renormalized by the subtraction of G(x, t) evaluated at 
x = t = 0, 



G R (x,t) = G[x,t) - (7(0,0) 



(27) 



This is equivalent to mass and coupling constant renor- 
malization. 
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Finally we evolve the field and two-point function ac- 
cording the equations (|2l|) and ( p2| ) on a discrete one- 
dimensional lattice. We employ a variety of lattices cor- 
responding to different physical volumes and lattice spac- 
ings and use a fourth order symplectic integrator to step 
in time. The range in the number of lattice points was be- 
tween 256 and 2048, with typical lattice spacing of 1/64 
to 1/128. Diagnostics performed at regular time inter- 
vals give a series of cinematic snapshots of the evolution 
and measure the components of the energy momentum 
tensor (T M „) as well as the separate components of the 
energy (Too), i-e. potential, gradient and kinetic energy. 

The results for the evolution of the quantum fields are 
shown in the next section, where they are also compared 
to the purely classical field bubbles of Section ||. 

IV. NUMERICAL RESULTS 

In Section |n| we have already discussed the paramount 
importance of Lorentz symmetry in constraining many 
of the properties of bubble wall propagation. These con- 
straints clearly also apply in the quantum theory at zero 
temperature. 

The first clear difference is that the field profile corre- 
sponding to the bubble is no longer purely classical. The 
quantum modes respond to our (initially) classical back- 
ground and partially screen it. The result is a dressed 
bubble, an object with a quasi-classical profile around 
which the vacuum is disturbed. This self-consistent ob- 
ject is Lorentz invariant and once formed propagates adi- 
abatically, i.e. without further particle creation. Sev- 
eral snapshots of the mean field and 2-point function in 
Fig. U illustrate this behavior, from initial formation to 
late propagation. 

The next important effect of quantum fluctuations is 
essentially static. Looking at Fig. || we also note that 
while the "height" of the dressed bubble (i.e. the true 
vacuum expectation value (vev) of the field) is initialized 
at the classical vev, it will not remain at that value dur- 
ing bubble growth. This is nothing remarkable, only a 
demonstration of the well known fact that the quantum 
effective potential differs from the classical potential, and 
their minima do not coincide. It might make more sense 
to initialize the field with the true vacuum expectation 
value given by the quantum effective potential at some 
loop order. However, given the naive initialization of the 
two-point function that we employ here, the field will 
still be perturbed away from this minimum inside the 
bubble when the two-point function changes. The ana- 
lyticly computated quantum effective potential treats the 
mean field as spatially homogeneous. Therefore it cannot 
be expected to describe accurately the dynamics of the 
bubble. We discuss this problem and its solution in more 
detail below and in Section V. 




X X 



FIG. 8. Snapshots of the expectation value for classical and 
quantum field (lower curves) and two-point function G (upper 
curve) at increasing time intervals. The initial configuration 
shown in the top left is almost identical in both cases. The 
evolution shows the quantum bubble relaxing into the mini- 
mum energy configuration preferred by the quantum effective 
potential while being screened by quantum fluctuations. In 
so doing, it leaves behind excitations inside the bubble and 
accelerates more rapidly. For presentation, G is shown shifted 
away from its actual value (equal to zero where <fi = 0.) 
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FIG. 9. The quantum Hartree approximation trajectories 
for initial bubbles of various sizes are shown as pluses. The 
lines are interpolations of the classical trajectories shown in 
Fig. [| The initial acceleration of the quantum bubble wall is 
manifestly more rapid, and it approaches its lightcone asymp- 
tote more quickly. 

At zero temperature, we observe both super- and sub- 
critical bubbles in the quantum case which grow and de- 
cay respectively as in the classical case. We track the 
position of the bubble wall in time. Despite the fact that 
the quantum evolution introduces fluctuations absent in 
the classical case, we see a bubble wall trajectory which 
is consistent with the classical picture. We observe ap- 
proximately hyperbolic trajectories which asymptotically 
approach the speed of light, as shown in Fig. |j| Fig. || 
shows some of the quantitative differences between the 
classical and the quantum evolution, including the ef- 
fect of the "roll" down to the minimum of the quantum 
effective potential. As the field vev changes inside the 
bubble, it releases potential energy beyond what would 
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happen from the movement of the interface-namely bub- 
ble growth-alone. Since this energy is converted into 
kinetic and gradient energy of the field and fluctuations, 
this roll can be held responsible for a more rapid accel- 
eration at the initial stage of growth. The energy of the 
true vacuum in the quantum case is also lower, so the rate 
of energy release by the bubble growth is larger. Hence, 
we see the quantum trajectories in Fig. ^| approach their 
asymptotic value (again unity) more quickly. This is the 
only significant departure from the particular hyperbolic 
trajectory predicted by a classical analysis. 

We also note that while the classical critical radius 
for growth was found to be R c = 0.35, the critical ra- 
dius in the quantum evolution is around R c = 0.45. In 
both cases, the critical radius for growth is a factor of 
4 to 5 smaller than the bounce radius Rb while rela- 
tively close (within a factor of 2) to the lower limit in 
(l+l)-dimensions for a wall with some thickness yT 1 (i.e. 
R c ~ 1/fi). In Section H we found /i" 1 = 0.24. 

It is worth considering whether the accelerated initial 
growth and the observed change in R c can be understood 
just from a classical analysis of the bubble using the effec- 
tive potential instead of the bare potential. Recall that 
the particular shape of the hyperbolic trajectory is deter- 
mined in the thin-wall analysis by the shape of the po- 
tential through the parameter Rb- For example, we can 
compute the effective potential numerically using homo- 
geneous fields and solving a gap equation along the lines 
of earlier works p6[ . From this, we can extract a quasi- 
classical potential and recalculate Rb and /i" 1 for this 
potential. We obtain Rb = 0.86 and /i -1 = 0.25. We 
fit the observed value of Rb = 1-3 from the trajectories 
in Fig. |^. While the classical thin-wall predictions using 
the quantum effective potential give the correct qualita- 
tive changes, they predict the wrong magnitude of the 
correction. We should also note that parametrically the 
thin wall approximation ought to be worse in the quan- 
tum case since now /iR ~ 3.44, less than half its classical 
value. 

This rough calculation harkens back to several efforts 
at applying the bounce analysis to quantum field models 
with, for example, symmetry breaking due to radiative 
corrections J27J. The problem in effect stems from the 
fact that the degrees of freedom which need to be traced 
over in the calculation of the effective potential cannot be 
properly integrated out since they participate in the bub- 
ble dynamics. We believe that a way out of this dilemma 
necessitates calculation of a self-consistent bounce, i.e. 
not only the shape of the mean field but also the full 
spectrum of interacting fluctuations in its background as 
prescribed by an effective action. While this is beyond 
the scope of the present work, we are pursuing it for fu- 
ture publication. Furthermore, with this information in 
hand, one could hopefully also bring to bear an analy- 
sis of the dynamical viscosity experienced by the moving 
wall at finite temperature, as suggested by other recent 



results 

In the meantime, the observed differences we have de- 
scribed in the shape and trajectory of the purely classical 
bubble and its dressed quantum counterpart are all that 
one may expect. At zero temperature, the qualitative 
characteristics of bubble propagation remain determined 
by the constraint of Lorentz invariance. Once finite tem- 
perature is considered, or alternately in the presence of 
a Lorentz-violating condensate, bubble propagation can 
change drastically, and the velocity of domain growth will 
in general not asymptote to that of light. 

Turning now to the last of Coleman's open questions, 
we consider the effect of colliding bubbles in the classi- 
cal and quantum Hartree approximations. Due to the 
nature of our numerical simulations, we enforce periodic 
boundary conditions on our one-dimensional lattice at 
the endpoints. When the bubble interface reaches the 
end of the lattice, it effectively meets its mirror image 
or, equivalently, a neighboring bubble wall. Whether the 
interfaces bounce off of each other or coalesce can then 
be observed. 

We note a surprising result: classically the bubble walls 
do not coalesce. The collision of the interfaces does tem- 
porarily excite wavelengths near the bubble "penetration 
depth" (the bubble thickness) at low amplitude, but this 
energy appears to be re-exchanged with the bubble inter- 
face. Furthermore, this is a special feature of the bounce 
profile which is not robust under deformations of the 
shape of the bubble: other initial bubbles do radiate and 
decay albeit slowly through collisions, see Fig. [U]. We 
also verified that the stability of the bounce disappears 
in (2+l)-dimensions where collisions rapidly destroy the 
spherical symmetry of the bubble. 

The (l-l-l)-dimensional classical bounce solution ex- 
hibits truly solitonic behavior in the sense that it is unal- 
tered by scattering off of another (solitonic) bubble. This 
effect, well known to exist in integrable nonlinear sys- 
tems [^o| , is unexpected in our model. In fact, Lohe has 
considered the linearized perturbations around the static 
soliton ( |fT|) of this model as one example in a parameter 
family of polynomial potentials (the sine-Gordon model 
appears as a limiting case); he has shown that there are 
completely reflected states in one direction from the soli- 
ton or bubble wall [|T|. This is just a physical conse- 
quence of the difference in the mass of elementary exci- 
tations about the two vacua: low-lying states in the false 
vacuum cannot pass through the interface because they 
have no energy counterpart on the other side. We may 
take this as suggestive evidence that the two bubble in- 
terfaces must repel under these circumstances, although 
it cannot guarantee the observed almost perfect reflec- 
tion. If the sign of the 7-term in the potential (Q) is 
flipped, the initial bubble will collapse and the interfaces 
will pass through each other, exploiting the cp — > —(f> 
symmetry of the model as they emerge on the other side 
of the collision region.. We verified this behavior nu- 
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merically. The classical scattering of our interfaces is 
thus quite analogous to the scattering of solitons in the 
sine-Gordon model where both types of behavior (perfect 
reflection and perfect transmission) obtain even though 
the vacua are identical |32fl . In that case, it is soliton- 
soliton or soliton-antisoliton solutions which display the 
two distinct possibilities. 

Quantum bubbles by contrast appear to dissipate en- 
ergy rather efficiently during collisions and especially 
during the propagation of the interface through the fluc- 
tuations created by the collision. Although at the level 
of our approximation the bubbles still rebound initially 
from the collision, they quickly lose kinetic energy. By 
the second or third collision, all semblance of the initial 
configuration is lost. In order to understand the quan- 
tum, or at least semi-classical scattering of the interfaces, 
one could perform a detailed study along the lines of 
Rcf . Q , adapted for the self-consistent quantum evolu- 
tion. 
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FIG. 10. Trajectories for bubble interfaces bouncing off of 
neighboring interfaces. Connected curves show classical tra- 
jectories for the bounce (bottom) and distortions thereof (up- 
per curves), while quantum trajectories are shown as data 
points. The classical bounce and the quantum data are plot- 
ted for lattice spacings of 1/64 and 1/128, but the trajecto- 
ries overlap entirely. The quantum trajectory loses coherence 
quickly relative even to the distorted classical bubbles. 

We present results in Fig. [To], where we plot the classi- 
cal and quantum trajectories at two different lattice spac- 
ings (same physical volume) in order to demonstrate that 
this effect is not an artifact of the lattice discretization. 
We also plot classical trajectories of distorted bubbles to 
show that they decay on a much longer time scale than 
the quantum bubbles. The trajectories are obtained by 
following one field value near the true vacuum; in the 
quantum and classical distorted cases, the apparent loss 
of periodicity indicates that fluctuations around the true 
vacuum no longer have a coherent structure. In Fig. ^ 
we plot a late time snapshot of the field evolution in the 
classical and quantum cases. 
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FIG. 11. Field expectation values for the quantum (up- 
per) and classical (lower) evolution after twelve collisions with 
neighboring bubbles (effected by means of periodic boundary 
conditions). The classical bubble still maintains its shape, 
while the quantum bubble appears to have dissipated its en- 
ergy into fluctuations around the true vacuum. 

V. DISCUSSION AND CONCLUSIONS 

In this paper we took the first steps towards showing 
how self-consistent quantum fluctuations can be incor- 
porated in the real time dynamics of bubble interfaces. 
We have shown that a classical bubble profile becomes 
dressed by quantum fluctuations, which in turn affect the 
rate of conversion of the false to true vacuum, i.e. the 
acceleration of the bubble velocity towards the speed of 
light. We have also found that the presence of quantum 
fluctuations promotes substantially more efficient trans- 
port of the bubble wall energy into particles at bubble 
collisions. 

There is one inelegant feature of our analysis, namely 
that the initial conditions for the quantum evolution were 
informed by the classical bounce calculation in absentia 
of self-consistent fluctuations. There should exist a so- 
lution to the coupled Dyson-Schwinger equations which 
is analogous to the bounce, the analytic continuation of 
which amounts to a dressed bubble which expands with- 
out radiation (particle creation). Using this configuration 
as a starting point for an initial value problem would 
eliminate the unwanted radiation that we observed re- 
sulting from the early time response of the vacuum to 
the nontrivial classical background. 

Starting from the 2PI-CTP formalism it is natural to 
generalize the one-loop computation of the false vacuum 
decay rate to a self-consistent problem, where fluctua- 
tions and the mean-field bubble profile are solved to- 
gether. The self-consistent nucleation rate and bubble 
evolution would be accurate to the same level of trun- 
cation as the effective action for any approximation, not 
limited to the Hartree example above. We believe that 
this intermediate step is necessary before more subtle 
features of the long time evolution-inevitably entangled 
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with the early time response-can be analyzed quantita- 
tively in detail. We have made progress towards this 
implementation, which we intend to report in a future 
work || . 

Knowledge of the spectrum of self-consistent fluctua- 
tions will render the inclusion of medium effects such as 
finite temperature more realistic. It also paves the way 
for other generalizations of the model and improvements 
upon the approximation we have employed. A finite tem- 
perature analysis in particular may be sensitive to higher 
loop order interactions which give rise to viscosity effects 
in linear response around equilibrium. Through analysis 
of the energy-momentum tensor (T^„) we can investi- 
gate the emergence of hydrodynamic interface propaga- 
tion from microscopic physics. We also intend to extend 
the model to include fermions. 

There are many more directions in which a similar 
analysis may be informative. From numerical solution 
of a subcritical bubble which decays, we can read the 
spectrum of asymptotic states. It may be possible that 
the same states, upon time reversal, would provide the 
initial conditions necessary to generate a bubble. Such 
an asymptotic configuration is generically a complicated, 
correlated many-particle state; however its overlap with 
a two-body state may hint at the probability of a bubble 
resonance in a scattering experiment. Stable dynamical 
solutions such as breathers have long been features of 
classical field theories. It is unknown whether they are 
stable upon inclusion of dynamical quantum effects. Fi- 
nally, there has been a great deal of excitement recently 
about extended objects such as branes and their dynam- 
ics, collisions etc. If there are fields confined to each 
brane in a homogeneous manner, then the effect of col- 
liding two branes can be cast into a problem very similar 
to the one explored here with the addition of traces over 
transverse degrees of freedom. 
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